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Abstract. We study the reduced dynamics of interacting spins, each coupled to 
its own bath of bosons. We derive the solution in analytic form in the white-noise 
limit and analyze the rich behaviors in diverse limits ranging from weak coupling 
and/or low temperature to strong coupling and/or high temperature. We also 
view the one spin as being coupled to a spin-boson environment and consider the 
regimes in which it is effectively nonlinear, and in which it can be regarded as a 
£^ ■ resonant bosonic environment. 

cr 

1. Introduction 

Comprehension of the phenomenon of decoherence in open quantum systems has 
always attracted much attention, in particular as a prerequisite to understand the 
transition from quantum to classical behavior. The dissipative two-state or spin- 
isj boson model has been thoroughly studied in wide regions of the parameter space with 

diverse methods and techniques since the 80's [UH]- 

In the last decade, the subject of decoherence has experienced renaissance 
following the growing interest in the field of quantum state manipulation and quantum 
computation [3J. Any noise source sensitively leads to a narrowing of the quantum 
coherence domain. This entails severe limitations for coupled qubits to perform logic 
quantum operations. For this reason, extensive understanding of the decoherence 
mechanisms is indispensable. 

In this work, we focus upon a model which is a generalization of the single spin- 
boson model to the case of two spins which mutually interact via an Ising-type coupling 
and are coupled to independent environments made up by bosons. The first analysis of 
this model relying on the influence functional method was given by Dube and Stamp 
|3]. They obtained results for the dynamics in analytic form in restricted regions 
of the parameter space by omitting certain classes of path contributions and bath 
correlations. Several other previous studies on the same or related models relied on 
the master equation and/or perturbative Redfield approach [SJ [7]. Besides the 
weak-coupling assumption, often the secular approximation 8 is made, which breaks 
down however when the spectrum becomes degenerate. 

| Author to whom any correspondence should be addressed. 
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The model allows, for instance, to study decoherence and relaxation of two 
coupled qubits [7] , or the influence of a bistable impurity on the qubit dynamics [§] ■ 
The latter may significantly degrade coherence in Josephson phase qubits [TU] . Other 
possible application is study of coherence effects in coupled molecular magnets [TT] . 

In earlier works, the model has been analyzed in the pure dephasing regime both 
by the Feynman- Vernon method [12] and the Lindblad approach [13] . Here we extend 
the work in Ref. [12j beyond the pure dephasing regime and include the full dynamics 
of the qubit. In particular, we are interested in the competition between decoherence 
and relaxation to the equilibrium state. Here we focus on the white-noise regime. We 
shall derive the exact solution for the reduced density matrix without restriction on 
the parameters of the model and analyze it in the coherent and incoherent domains 
and in the crossover regions in between. 

The model and relevant quantities of the reduced dynamics are introduced in 
section 2. Section 3 deals with the exact formal solution for the reduced dynamics. In 
section 4, the path sum is carried out in the white-noise domain without any further 
approximation, and analytic expressions for the relevant expectation values in Laplace 
space are presented. After an overview of the qualitative features of the dynamics in 
section 5, we present in section 6 explicit expressions for decoherence and relaxation in 
the various parameter regimes ranging from low temperture and/or weak coupling to 
high temperature and/or strong coupling. Finally, we study in section 7 the influence 
of a nonlinear spin-boson environment on the second spin in the various limits. We 
demonstrate that it behaves in the weak-coupling limit as a bosonic (linear) bath with 
a resonant spectral structure. 

2. Model 

We consider two two-state systems which are coupled to each other via an Ising-type 
coupling and to independent bosonic environments. In pseudospin representation, we 
choose the generalized spin-boson Hamiltonian (we use units where h = fc B = 1) 

H = — y Ux ~ ~T Tx ~ 2 azTz ~ 2 azXl ~ 2 TzX2 + ^ H w o & L 6 o ■ W 

C=l,2 a 

In the basis formed by the localized eigen states \R> and \L> of a z and t z , respectively, 
Ai and A2 represent the tunneling couplings between the localized states, and the 
coupling term — \va z r z acts as a mutual bias energy of strength v. The collective bath 

modes X^(t) = J2 a c c ai^c «(*) + ^ a Wl (C = 1, 2) represent fluctuating bias forces. 
The Hamiltonian is very rich in content and may model diverse physical situations. It 
may describe two coupled qubits or a qubit a in contact with a complex environment 
formed by a bistable dissipative impurity r. Other possible realizations are coupled 
molecular magnets of which the low-energy states can be viewed as a spin [llj . 

For the model |T]), all effects of the environments are captured by the power 
spectrum of the collective bath modes 

S cc (cu) = ± J™dt (X c (t)X cl (0) +X c (O)X c (t)) = <y u ,7rG c («)coth(^) (2) 

with the spectral density of the coupling [TJ [5] 

G c (u;) = ^2<? Ca 6(w-w<. a ) = 2tf c we-W/'"«, C = 1, 2 . (3) 
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Here, the second form represents the Ohmic case with a high-frequency cut-off ui c . 
Alternatively, one may choose that the two spins are coupled to a common bath [TJ. 
Here we study the effects of independent environments. This case is realistic in most 
physical systems of actual interest. 

The density matrix of a single spin has four matrix elements, the two populations 
that we shall label as RR = 1 and LL = 3, and the two coherences with labels LR = 2 
and RL = 4. The two-spin density matrix has 16 matrix elements p nm {t). We 
choose for convenience that the first (second) index refers to the states n = 1 , • • • , 4 
(m = 1, •■•,4) of the cr-spin (r-spin). The matrix elements can be expressed in 
terms of expectation values of 15 operators, (<7j <g> l) t — (<Ji)t, (1 <8> Ti)t — (n)t and 
((Tj ®Tj) t — (o~iTj) t (i — 1, 2, 3 and j = 1, 2, 3). The 4 pure populations may then be 
written as 

Px,x(t) = [1 + (*z)t + (T z )t + (<r s T z ) t ]/4 , 

Pi, 3 (*) = [1 + (O* - <7i>t " <<V a >t]/4, 

(4) 
P3,i(*) = [1- <Ot + <r,) t - (o-,r,) t ]/4, 

P 3 , 3 W - [1 - K)* - (r,)* + (%)t]/4. 
Corresponding expressions hold for the 4 pure coherences and the 8 hybrid states. For 
instance, we have 

02,4 (*) = {( a x T x)t + ( a y T y)t + H^x T y)t - ' l { a y T x)t ]/4 ■ (5) 

Here we are predominantly interested in the populations. Throughout we will choose 
that the reduced system starts out from the initial state p l 1 (t = 0) = 1 while the heat 
reservoirs are in thermal equilibrium at temperature T. 

In the absence of the environment, the Hamiltonian H — Hq can be easily 
transformed into diagonal form 

O A 

H = — - (a z ® 1) - - (1 ® t z ) . (6) 

The eigenfrequencies are 

q = i(n+ +n_), 5 = l(n+ -n_), (7) 



with n± = V(Ai±A 2 ) 2 +w 2 , (8) 

and they obey the Vieta relations 

n 2 + 5 2 = Aj + Al+v 2 , n 2 s 2 = A 2 A 2 , 

(9) 

n 2 _ + n 2 + = 2(a 2 + a 2 + i; 2 ), n 2 + n 2 _ = (A 2 + A 2 + v 2 ) 2 - AA 2 A 2 . 

The Liouville equations Wj(t) — —i[H, Wj(t)] (j = 1, - - • , 15), where the set 
{Wj(t)} represents the above 15 operators, yield 15 coupled equations. These are 
conveniently solved in Laplace space. For instance, we get 

A(z, 2 + A 2 + A 2 ) A( V 2 + A 2 + A 2 ) 



(A 2 + n 2 )(x 2 + s 2 ) ' " z v " (a 2 + n 2 )(\ 2 + s 2 ) ' 

, ms (v 2 + X 2 )(v 2 +A 2 + A 2 + X 2 ) (10) 

V* T *Wl - A(A 2 + f> 2 + )(A 2 + f> 2 _) 
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Here we are interested in the evolution of the two-spin system without restricting 
ourselves to weak damping. Therefore, we refrain from employing the perturbative 
Redfield approach. Rather we calculate the reduced dynamics with use of the 
Feynman- Vernon influence functional method. We show that the solution is available 
in analytic form in the white noise limit for general parameters A 1; A 2 , v and T. 

3. Formal solution for the reduced density matrix 

Within the Feynman- Vernon method, the exact formal expression for the RDM of the 
two-spin system is the quadruple path integral 

P ai<T ,. T(T ,{t) = IvoVo'VtVt' A[a}A*[<T'}A[T]A*[T'}B[a,o-';T,T'}f[a,o-';T,T'} (11) 

with appropriately chosen boundary values for the spin paths. Here, each of the 
paths cr(t'), cr'(t'),r(t'), r'{t') starts out from the localized state \R> at time zero. 
They end up at time t in the states |c f >, |<7f>, |T f >, and |r f '>, respectively, where 
er f , fjf , r f , Tf £ R, L. The functional A[cr] is the amplitude for the free spin a to follow 
the path a(t'), the functional B[a, a'; t, t'} represents the coupling of the two spins (see 
below), and the functional J r [a,a']T,T l ] introduces the environmental influences. 
For uncorrelated baths, we have T\o,o'\t,t'\ = .Fi[£i,?7i] ^2^2! %]> where 

Jo Jo L J 

Here we have introduced symmetric and antisymmetric spin paths, 

ei(0 = \W)-a\t')], Vl (t') = \W)+a\t')], 

2 2 (13) 

^ 2 (*') = ^W)-r'{n, V2(t') = -[T(t')+r'(t')]. 

The correlator Q^(t) = Q'At) + iQ'i{t) is the second integral of the force auto- 
correlation function (XAt)XA0))p (see eq. ©). In the Ohmic scaling limit, we have 



Qc(t) = 2K C 



In | — - I + In sinh I — ^ 



iirK c sgn(t) , C = 1, 2 . (14) 



Here, Kq is the usual dimensionless Ohmic coupling strength for the spin £, and 
/3 = 1/T is the inverse temperature. 

To handle the quadruple path integral (fTI]) , we follow the procedure for the single 
spin-boson problem [H[2] and write it as an integral over two paths, one for each spin. 
Each such path visits the diagonal "sojourn" states and the off-diagonal "blip" states 
of the respective spin. A path which starts and ends in a sojourn state must contain 
an even number of transitions with amplitude =FiA^/2 for each flip of spin £. The 
flips occur at times i ■ for spin 1 and at times s for spin 2. Upon labeling the sojourn 
and blip states with charges n 1 j, ^ • (for spin 1) and r] 2 •, £ 2 ,j (f° r s P m 2), each with 
values ±1, the paths with 2n x and 2n 2 transitions, respectively, may be written as 

«i 
rit H \t') = J2vi,j W t 2j ) 6{t' - t 2j+l )} , (15) 

3=0 
"l 

el" 1 V) = E^j ^' ^-^ - e{t ' - i2 ^] > ( 16 ) 
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n 2 

3=0 

n 2 

&>\t') = £e 2j W S2j _ x ) 6(t' - s 2j )} . (18) 

Upon introducing the notation Q(-,j,k — Q( {tj—tk), we may write the bath correlations 
between the blip pair {j, k} of spin £ in the compact form 

h-(;j,k = Q(;2j,2k-1 + Qc;2j'-1,2& ~~ Q(;2j,2k ~ Q C,2j -l,2k-l ■ (19) 

With this, the influence functional for the paths (fT3"]) reads 



_ ' l c j'-i 

^c"' ' : ex P [ - £ <9c;2j,2i-il exp [ - J2 £ £cA* A to,* 

j=2 fe=l 



x exp 



3 = 1 


exp ^ 


rif — 1 
i7T^ c E Cc.fc+l^C* 



(20) 



fc=0 

Here the first and second term represent the intrablip and interblip correlations, 
respectively. The phase term is specific to the Ohmic scaling limit and represents 
correlations of the sojourns with their subsequent blips. 

The sum over all paths now means (i) to sum over all possible intermediate sojourn 
and blip states of the two spins the paths with a given number of transitions can visit, 
(ii) to integrate over the (for each spin) time-ordered jumps of these paths, and (iii) 
to sum over the possible number of transitions the two spins can take, 

£•••-£ E E / d *2% / d Vi- di i 



all paths n 1 =0{^ 1 _ 3 =±l}{5 1J =±l} 



(21) 



X E E E / dS 2n 2 / ds 2n2 _ 1 --- d Sl ---. 

™2=0{%, 3 =±i}U 2>3 =±i} J * 1 J *. J *. 

For coherences, the number of transitions in the respective spin path is odd. 

Next, we rewrite the double path sum (|21[) with time-ordering for each spin in 
terms of a single path over the 16 possible states with time ordering of all the flip times 
{£-•}) {sj}. In this representation, the system lingers for some period in a particular 
state of the double-spin RDM and then it flips with amplitude ±i A x /2 or ±i A 2 /2 
to another state. For the longitudinal spin-spin coupling ~^v<t z t z in ((TJ, the coupling 
factor B in eq. ((lip is unity when the system dwells on one of the four pure soujourn 
or one of the four pure blip states, and it is sensitive to the coupling when it stays, 
say for a period u -, in one of the eight hybrid states. In detail, we have 

B[j] - e +iv "' , for j e (1,2), (2,1), (4,3), (3,4), 

(22) 
B[j] - e""^, for j e (1,4), (4,1), (2,3), (3,2). 

Combination of the above expressions yields the exact formal solution for the 
dynamics of the RDM of the two-spin system in the Ohmic scaling limit. Evidently, 
because of the nonconvolutive form of the bath correlations in the influence functional 
(|20p. the path sum can not be performed in analytic form. Alternatively, one may 
recast the exact formal series expression for the populations in the form of generalized 
master equations in which the kernels, by definition, are the irreducible components 
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of path segments with diagonal initial and final states [2]. In the general case, the 
kernels are given by an infinite series in A x and A 2 , with the time integrals in each 
summand being again in nonconvolutive form. Additional difficulties in performing 
the path sum (|2ip arises from the spin-spin coupling 



4. Exact solution in analytic form in the white-noise limit 

In the white-noise limit of eq. @, S+ (u> -C 1//?) = 2$*, where 

i? c = 2ttK c T (23) 

is a scaled thermal energy, the bath correlation function takes the fornj§| 

Q c (i) = 2tf c [ln(^_) +T C ] +0 c \t\ + inK c sgn(t) . (24) 

This expression emerges directly from eq. (|14p in the high temperature or long-time 
limit t/(3 ^> 1. The first term in eq. (l2"4l> leads to an adiabatic (Franck-Condon-type) 
renormalization factor made up by modes in the frequency range 2t:T < lo < lo c . 
It is natural to assimilate this term, together with the phase term, into an effective 
temperature-dependent tunneling matrix element, 

A 2 = A 2 T cos(nK c ) , with A 2 , T = (27rT/A c , r ) 2X c A^ r e- 2 ^ T < , (25) 

where A f r c = A^/oj c c is the standard renormalized tunneling matrix element. 

All dynamical effects of the environmental coupling are captured by the second 
term AM in eq. (|2"4]l . From this we see that the weight of the thermal energy relative 
to the systems energies A, and v is assessed by the scaled thermal energy #>. 

Based on experience from the single spin-boson system we should expect that the 
form (|24p is a remarkably good approximation in the parameter range [2] 

K c < 0.3 , and Q± < T < lo c , (26) 

and this is corroborated indeed by our study. For a single unbiased spin with Ohmic 
damping K <C 1, the coherent-incoherent "phase" -transition is at temperature T = T* 
with T* = A T /(nK) [2]. Therefore we should expect, that, for K^ < I, the white 
noise form (|24p is valid not only in the incoherent regime but also in a sizeable domain 
of the coherent regime. This shall be confirmed subsequently. 

For the form (|24p of QAt), the interblip correlations cancel out exactly in eq. (fiT))) , 
A,- , i, —Q. As a result, each term of the infinite series for the RDM becomes a convo- 
lution. This makes the path sum accomplishable. We shall now exemplify, by taking 
(<r z )t and (a z T z ) t as examples, that the path sum for the Laplace transform of the 
RDM can be carried out exactly in analytic form. This is achieved by first calculating 
the kernels and then summing up the respective geometrical series of these objects. 

4--1- The expectations (o~ z )t and {r z ) t 

By definition, the kernels represent irreducible path segments which interpolate 
between pure sojourn states. Irreducibility means that these segments can not be 
separated into uncorrelated pieces without removing bath correlations. The analysis 
gives that every contribution to the kernel /C(A) of (a z (\)) displays initially and finally 
a transition of spin a with any number of even hops of spin r at intermediate times, 

§ The term T^ accounts for deviation of the actual high frequency behavior of G^(u>) from the 
exponential cut-off form in eq. J3) [2] . 
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as shown in the diagrams of Fig. [TJ There are no other contributions. For all other 
irreducible diagrams, one may think of, e.g., where either the first or the last flip, or 
both of them, are flips of spin r, the respective contributions from the two different 
final states of spin r cancel each other. 






Ai A 2 



Figure 1. Sketch of the irreducible kernels (with amputated legs) /Ci, (A) 
(a), /Ci i i(A) (b), and ICi^(X) (c). The solid and dashed curve represent the 
correlations of bath 1 and of bath 2, respectively. 

It is convenient to write /C(A) = J^Lo^i n(A); where IC ln (X) is of order A 2 A 2 ,™. 
The sum over all spin states of order (1, n) yields the expressions 



^i.i(A) 



Af T / dr : e 



- (Ai '' ir i cos(7rif 1 +vt 1 ), 

A 2 T A 2 T / dr 1 dT 2 dr 3 e-( A+ ''i)^ + ^+^)e-^^ 
Jo 
1 

COS^J^ + VT X — VT 3 ) ] C0S(7T K 2 ) 



x - [cos(ttK 1 



VT-, 



VT. 



(27) 
(28) 
(29) 



Here we have taken into account that, for the white- noise form (|24|) . a correlation of 
bath £ stretching over an interval between neighboring hops effectively leads to a shift 
of the Laplace variable in the respective time integral, A — ► A + i), . It is convenient 
to split the kernels into the contributions which are even and odd in the coupling v, 
/C(A) = /C' + )(A) +fc(~\\). The resulting expressions may be written as 

A + ^ 



Ki,„(A) = K M (A) -A| / dr^e-^+^e-^+^cos^) 



/C(+)(A) = -A 2 



(X + ti^+v 2 



\2 



(30) 



(A + ^)(A + ^ + 2 ) ( A + tf i)2 +V 2 + -^-Aj 



/C(-)(A) = tan^iv^A 2 



A + ^ 



A 2 

-*2 



(31) 



A + tfi+A, (A + ^) 2 + « 2 + i^%lAl 



Paths which visit a pure sojourn state at intermediate times yield reducible contri- 
butions. Taking into account all possibilities of such visits yields a geometrical series in 
the kernel /C*- + '(A), while IC^'(X) occurs only once as initial irreducible contribution. 
Thus we get for (cr z (A)) the concise form 



= I I + /C<-)(A)/A _ 1 A^A) 
^ * W ' A 1 - /C(+)(A)/A A D t (A) 



(32) 
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Algebraic manipulation gives (o~ z (\)) in the form of a simple fraction with denominator 
and numerator in the form of polynomials of degree four. We get 

JV X (A) = A 4 + (30 X +0 2 )A 3 + (v 2 + A 2 + 3-d 2 + 2$ 1 $ 2 )\ 2 

+ {v 2 d 1 + Al^ + d\ + v 2 -d 2 + tf?tf 2 )A + vA 1 t3.n(nK 1 )(d 1 + # 2 + A) , 

D X (A) = A 4 + (S^+^X 3 + (fi 2 + 5 2 + 3tf? + 2tf 1 tf 2 )A 2 

+ (v 2 ti 1 + 2A 2 d l +A 2 i) l +i)l+v 2 ti 2 + A 2 d 2 +d 2 d 2 )\ (34) 

+ A? -&\ + A 2 i? x # 2 + n 2 s 2 . 

Here we have introduced the eigenfrequencies Q, S of the undamped coupled two-spin 
system. The bar denotes adiabatic renormalization in eq. according to eq. (|25|) . 
The pole at A = in eq. ([32| yields the equilibrium value 

/CC-)(Q) tt MgjjrA) f o^ 

^./.q /c<+)(0) 2T A2 + tf 1 (i? 1 + i? 2 ) ■ [ ' 

Hence (c z ) is negligibly small for A 2 7^ 0, while, as A 2 — * 0, it takes the proper 
(white noise) equilibrium value v/(2T) of the single biased spin boson system. 

The dynamical poles A ■ (J = 1, • • • , 4) are given by a quartic equation with real 
coefficients. Upon collecting the various pole contributions, we get in the time domain 

4 

(a z ) t = Y. A ^ + (Oeq- ( 36 ) 

1=1 
Evidently, the expectation (r z (A)) takes similar form, 

{TAX)) ~ XDzJXJ- (37) 

Here, the polynomials iV 2 (A) and Z?2(A) follow from the expressions (f3"3"| and (j3"4"|) by 
interchange of the indices 1 and 2. 

(b) 




Ai A 2 



Figure 2. Sketch of the irreducible kernels ^^(A) (left) and ^ b) (A) 
(right). The intervals are dressed by self-energy contributions of spin a 
(black circle) and spin r (black square), as sketched in Fig. [3] 



A] A! A! A! A! A! 

• ~ X / ~ \ ' ~ \ 

-J l_ + —I u-L 1_ + .. . 

A 2 A 2 A 2 A 2 A 2 A 2 



Figure 3. Self-energy terms due to spin a (circle) and spin r (square) 
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4-2. The expectation (cr z T z ) t 

Consider first contributions to the kernel of (cr z r z (A)) in which the first and last flip 
are made by spin er, as sketched in Fig. O In the intervals of the bare diagrams, either 
spin a or spin r or both stay offdiagonal. Every interval in which spin £ dwells in 
a sojourn state is dressed by selfenergy contributions schematically given for spin a 
(circle) and spin r (square) in Fig. [3] Diagram[2](a) yields 

4 a) (A) = K x (A) V - A 2 \ &T X dr 2 e ~Wi>i e -Wi+»,)r> cos^) 

n=0 V Jo 

= - (A + ^ + i? 2 )[l- ^-^-tan( 7 ri<: i )]a 1 (A) ) 

where the functions ay (C — lj 2) are given by 
A + $ A 2 

MA) = ^^ (A + ^ + ^ + ^^Af ' (38) 



« 2 (A) = , , M 2 , a — — . 2 A+ , T „ ■ (39) 



A + 1? 2 A 2 

A + tf 1+ tf 2 (A + ^ 2 ) 2 +, 2 + ^±^ A? 

The nested diagram in Fig.[2](b) produces the additional factor a 2 (A)a 1 (A), 

4 b) (A) = 4 a) (A)a 2 (A) ai (A). (40) 

Higher-order nested diagrams of the type sketched in Fig. [2] class into a geometrical 
series. All these terms are readily summed up to the contribution 

A X (X) = 4 a) (A) 1 ^tt • (41) 

1 - a 2 (A)a 1 (A) 

Similarly, we find that diagrams (a) and (b) in Fig. @] yield the expressions 

B ( x a \x) = (\ + 4 1 + # 2 )[i--J^toi{irK 1 )]a 1 (\)a a (\), (42) 

B[>\\) = B< a >(A) ai (A)a 2 (A). (43) 

With all higher-order nested diagrams of this type added, one finds again a geometrical 
series in Q! 1 (A) a 2 (A), 

B X (X) = B[ a \\) - j—j- . (44) 

1 - Q 1 (A)a 2 (A) 

Clearly, we must add those contributions resulting from the terms A X (X) and B X (X) 
by interchange of the spins a and r. The analysis shows that there are no other con- 
tributions. Again, we split the kernel C(A) into the parts which are even and odd in 

the coupling v. We readily get for C^\X) = £ c=1)2 [ -^ ( A ) + B^\\) } the forms 

C^(X) = -(X + x + * a ) «i(A) + a 2 (A) - 2a (A)a 2 (A) 

1 - a^Ajo^A) 

, , ^-^ u av(A) — a-, (A)a 9 (A) 

C<"> A = (A + tf 1+ tf 2 ) ^ ^-^tan(7rX c )^ , , ■ (46) 

r— ^ A -I- w A s 1 — a 1 (A)a,(A) 



C=i,2^"C 



These expression represent the entity of irreducible path segments. Next, we observe 
that the sum of two-spin paths with any number of interim visits of pure sojourn 
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states yields a geometrical series of these objects. In the part which is odd in v, the 
first irreducible path section is again described by the kernel C^(\). Thus we get 

1 1 + C(~)(A)/A 1 N(X) 
{a * T * {X)) = A 1 - C(+)(A)/A = A D(X) ■ ( } 

The second form is a simple fraction with polynomials N(\) and -D(A) of sixth order, 
N(X) = w 3 (i? 1 + l ? 2 + A)[tan(7ris: i )A2 + tan(7rir 2 )A2] 

+ t;tan( 7 r J fC 1 )A2(i? 2 +A)[A2 - A| + (0 2 + A) (0 X + 2 + A) ] 

+ «tan(7riY 2 )A 2 (tf 1 +A)[A 2 - A? + (0 X + A) (0 X + 2 + A) ] 

+ Aw 4 (tf 1 +tf 2 + A) + Aw 2 [A2(^ 1 +A) + A2(^ 2 + A)] (48) 

+ Xv 2 (d 1 + tf 2 + A)[tf 2 + tf 2 . + 2A(0 X + tf 2 + A) ] 

+ A (0j + A)(0 2 + A)[ A? (0 X + A) + A 2 (0 2 + A) ] 

+ A(tf 1 +A) 2 (tf 2 + A) 2 (tf 1 +tf 2 + A), 

D(A) = A 6 + 3(?$»i+7? 2 )A 5 + (^ 2 h + fi 2 _+3^ 2 + 8^n? 2 + 3^ 2 )A 4 

+ (ft + t? 2 )(2^ + 2Q 2 _ +tf 2 + 6i?i0 2 + 01) A 3 

+ [O^O 2 . + 0x02(20! + tf 2 )(0! + 2tf 2 ) + v 2 {M\ + 40!0 2 + 30|) 

+ A 2 (2^ + i? 2 )(0! + 30 2 ) + A 2 (20 2 + X )(0 2 + 30 x )] A 2 (49) 

+ (0! + ^[ft 2 . ft 2 . + ^2^| + (t> 2 + £2)02 + (v 2 + £2)02 

+ 3^^ 2 (A 2 + A 2 )]A + W 2 (^ 1 + ^ 2 )(A 2 ^ 1 + A 2 ^ 2 ) 

+ 0!0 2 [ (A 2 - A 2 ) 2 + (A 2 ^! + A 2 ^)^ + 2 ) ] . 

(a) 





Figure 4. Sketch of the irreducible kernels Bj (A) (left) and B^ (A) (right). 
Again, the intervals are dressed as sketched in Fig. [3] 

The odd powers of the pole function can be removed with a shift. Putting 
D[ A = x — (i} 1 + 19 2 )/2] = D(x), we obtain a polynomial with even powers, 

d(x) = x 6 + [n 2 + + n 2 _- Ud 2 + d 2 2 ) -Wr- $ 2 ) 2 }x i 



(50) 



+ [ti 2 + n 2 _ + (A 2 + A 2 - 2v 2 )ti 1 d 2 - A 2 2 ~ A 2 2 
+ U^i ~ ^ 2 ) 2 (3^ + 20 x 2 + 30 2 ) ] x 2 

^- 4 [v 2 (^+^) + (a 2 -a 2 )(^ 2 -^) + m-^)(^-^)?- 

Thus we have in the time domain 

6 

(<r,T s ) t = Y. B * cKt + K r *>e q = (51) 
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where the {A •} are the zeros of -D(A) = 0, and the equilibrium value is 

e(->(o) _ v_ 



The expressions (|32[) - (|37j) and (j47|) - ([50)1 are the main results of this work. 
They represent the exact analytical solutions for (a z (X)}, (r z (A)), and (ct z t z (X)), in 
the white-noise limit for general coupling v and general effective reservoir couplings 
, d 1 and d 2 - Except for use of the form (j!M|) , no other approximation has been made. 
We remark that for all other initial and final states of the RDM we would find the 
same pole functions (j34|) and (|49[) . Only the numerator function would be different. 



5. Qualitative features 

The behaviors of the four dynamical poles of (o~ z ) t and the six dynamical poles of 
(o~ z T z ) t , and the respective amplitudes are quite multifarious. In this section, we 
sketch the characteristics for the symmetric system, Ai = A2 = A and $ 1 = § 2 — &• 

5.1. (a z ) t : 

In the coupling range v < v CI — A/V2, there are three crossover temperatures, denoted 
by #5, #£, and d* 2 (see Figs. E and E>). 

In the regime d < $\ the dynamics is coherent and described by a superposition 
of two damped oscillations. For d < $g, the oscillations have different frequency and 
the same damping rate, and the amplitudes are comparable in magnitude. On the 
other hand, in the range $q < $ < #*, they have the same frequency, but different 
decrement and the amplitude belonging to the larger decrement is negligibly small. 

In the temperature regime 1? > i?J, the dynamics is incoherent. In the regime 
$1 < •& < $2, the four poles are real, and the two smallest rates have largest amplitudes 
and dominate the relaxation process. 

In the so-called Kondo regime § > $|, the dominant pole is real and approaches 
— A I'd, and its residuum goes to 1 — (c 2 ) cq , as temperature is increased. The other 
real pole takes the value — 2$, while its residuum drops to zero. There is also a damped 
oscillation of which the frequency and rate approach asymptotically v and #, but the 
amplitude becomes negligibly small. The phenomenon that, in the Kondo regime 
for K < J) incoherent relaxation slows down with increasing temperature, is already 
well-known in the single spin-boson problem [2], 

Fig. [7] shows the transition from coherent to incoherent dynamics as $ is raised. 
One can also see that at high -d the effective damping decreases with increasing -d. 

For v > v CT , there is only one crossover. It separates the regime with two complex 
conjugate poles from the regime with one pair of complex conjugate poles and two 
real poles. Above $*, the relaxation is governed by the Kondo pole. 

5.2. (<r z T z ) t : 

For a symmetric system, two poles are cancelled. The reduced pole equation reads 

A 4 + MX 3 + (Cl 2 + + 5tf 2 )A 2 + 2tf(n 2 + + tf 2 )A + M 2 A 2 = . (53) 

The characteristic behavior of the poles as function of the scaled temperature $ is 
shown in Fig. [8] At low $, all poles contribute to the dynamics. In the Kondo 
regime i9 ^ 3A, the leading pole behaves as — 2A 2 /$, and the amplitudes of the other 
contributions are negligibly small. 
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Figure 5. (cr z ): Crossover temperatures as functions of v, A = 1. 
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Figure 6. (cr z ): Real (a) and imaginary (b) part of A ., v = 0.5, A = 1. 
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Figure 7. (u z ): Plots of (cr z ) in the Fourier regime (a) and time regime (b), 
v = 0.5, A = 1, K = 0.05. Red (small-dashed) ■& = 0.2, green (dashed) ■d = 1.2, 
blue (medium-dashed) i9 = 2.1, grey (long-dashed) i? = 2.8. 



6. Dynamics in the various parameter regimes for differing spins 

6.1. Low temperature behavior 

Below the first crossover temperature, $7 ± ^T < T *, the real parts of all poles are 
varying linearly with T. In detail, the poles and amplitudes are as follows. 
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Figure 8. (cf z t z ): Real and imaginary part of A , v = 0.5, A = 1. 



6.1.1. (o~ z )t'- There is a superposition of two damped oscillations with frequencies 

n 2 + A 2 -2P 



^1,2 
^3,4 



= ±i<S - 
= ±iO 



2{n 2 -P) 

2n 2 -8 2 -A 2 

2{n 2 -P) 



h 



O 2 - A 2 

"' 2 ^ 9 



01 



2(ft 2 -5 2 ) 2 

Aj-5 2 
2 (ft 2 - <5 2 ) 



tfo 



0(tf 2 ) , 
- 0{ti 2 ) 



and amplitudes 



^ = i%^i! + ow , 



^3,4 — - 



1 n 2 - « 2 - A 2 



ft 2 -<5 2 



(54) 



+ 0(t?) . (55) 



There is a superposition of two damped oscillations and two 



6.1.2. (a z T z )i 

contributions describing incoherent relaxation towards the equilibrium value (cr z T. 



A 



1,2 ~~ 
Ak = 



±ifi 



^ + tf 2 



v 3,4 



±m 



n 2 - A? n ft 2 - A 2 

-^i - 



A 2 -J 5 



o 2 - s 2 Ul " n 2 - s 2 "° 2 ' Xe n 2 - s 2 Ul n 2 - s 2 " 2 



^i 



Ai ~^ 



2/cq' 



(56) 



where terms of order ©(i? 2 ) are disregarded. In leading order, the amplitudes are 



^-iV-m 



B 3A = -( I 



B K 



v 2 S 2 



(n 2 -s 2 ) 2 n 2 



V 

4 l - W 

B R = 



2o2 



v 2 n 



{Q. 2 -8 2 f 



n 2 - s 2 



(57) 



(58) 



6. 2. The regimes of large coupling and/or high temperature 

When the coupling v and/or the scaled temperatures $ 1 2 are large compared to the 
other frequencies, the amplitudes of three (five) pole contributions to (cr z )t ((c z T z )t) 
are negligibly small, and among the dynamical poles only the real pole with smallest 
modulus is relevant. Hence the two spins essentially behave as a single spin which 
relaxes incoherently to the equilibrium state according to 



{Vzh 



-7„t 



z I eq ' 



i (J z T z)t = (1 - {^ Z T- Z ) eq )e 



'ctt 



(59) 
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6.2.1. (cr z ) t : The relaxation rate is found from -D 1 (A) = with the form (|M|) as 

= _&j_ A| + tf 1 (* 1 -H> 2 ) 
7ff v 2 + tf 2 ^ 1 + ^ 2 ■ W 

This reduces in the parameter regime , d 1 3> v, A 2 to 

la = A?M « T 2 ^- 1 . (61) 

In this regime, (o~ z ) t is independent of the coupling v and hence independent of the 
dynamics of the r-spin. The temperature dependence j a ex T 2K ^~~ l distinguishes the 
so-called Kondo regime, in which, for K 1 < i the relaxation dynamics slows down as 
temperature is increased. 

On the other hand, when v 3> A 1 2 , the two spins are locked together [4], and the 
effective tunneling matrix element is 6 = AjA^v, as follows from J2J with ((8]). We 
then get from eq. (f60|) 



s- 



2 



la = 



d 1+ d 2 



1 



*l(01+*») 



A 2 



(62) 



This yields the limiting expressions 



la = ^K (A 2 «^ lj2 ), and 7cr = _i-^_ (A 2 »^, 2 ). (63) 

The former is the relaxation rate of the biased single spin-boson system at low i? x . 
The latter describes Kondo-like joint relaxation of the locked spins. 

6.2.2. (a z T z ) t : In the incoherent regime, the relaxation rate r y aT of the effective 
single spin receives rate contributions from both the a- and the r-spin as if these were 
independent biased spins in contact with their own heat reservoir. 

In the large-coupling limit, v 3> A 1 2 , $ l 2 , the relaxation rate is found as 



A 2 A 2 

7ctt = -M + T-fV ( 64 ) 



The individual contributions are single-spin rates in the large-bias regime. In the 
high temperature limit $ l 2 3> v, A x 2 , on the other hand, both rate contributions are 
Kondo-like, 

A 2 A 2 

■<- " §f + t. ■ (65) 

Consider next the regime d x S> A l5 A 2 , v, in which spin a behaves Kondo-like, 
as in eq. (|6ip . Hence the dynamics of the a-spin is slow compared to that of the 
r-spin. Thus we should expect that (a z T z ) t approaches the dynamics of the biased 
single spin-boson case as $ 1 is increased. Taking into account terms of linear order in 
7 CT in the pole equation, the expression (|47| with (|48l) and (|49|) assumes the form 

A 2 + 2t9 2 A + v 2 + i? 2 + TrK 2 vA 2 /X 

{a * TA >> ~ A^ + 2{d 2+la )\ 2 + (A 2 +t, 2 +tf 2 +3 7ff tf 2 )A + A 2 tf 2 + 7 > 2 +tf 2 ) ■ 

Indeed, in the limit 7o . — ► 0, this form reduces just to the analytic expression for 
(r 2 (A)) of the biased single spin-boson system in the white- noise limit [2]. 
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Let us now view spin r with reservoir 2 as an environment for spin a . This complex 
environment is in general non-Gaussian and non-Markovian [14j . Recently, the same 
model has been studied numerically using a Markovian master equation approach [15J. 
To proceed, we first note that in the absence of bath 1, i^ =0, the pole equation 
N^X) = is still of fourth order. There is no reduction in the general case. In Fig. [9] 
we show plots of the four poles as functions of i? 2 f° r a particular set of parameters. 
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Figure 9. (<r z ) with spin-boson environment: Real (a) and imaginary (b) part 
of A . At low i? 2 , {u z ) t is a superposition of two damped oscillations. At high 
i9 2 , there is one damped oscillation and one relevant relaxation contribution. The 
parameters are v = 0.8, Ai = 1, A2 = 1.5. 



7.1. High temperature limit 

Simplification occurs, however, when i? 2 is very large compared to the other 
frequencies. In this regime, the kernel (|30|) reduces to the form 



JC£>(A) = A- A? 



l + ~ 



It 



A 2 + «H' A A 2 + v 2 +-y T \J (66) 

where j T = A 2 /-^ is the relaxation rate of spin r in the Kondo regime. With this 
high-temperature expression for the kernel, the quantity (a z (X)) is found to read 

1 A 2 + 7r A + i; 2 



MA)) = 



(67) 



A + 4^ (A) A3 + ^ + ( A ? + ^) A + A ^r 

This expression describes the dynamics of spin a coupled to a spin-boson environment, 
where the latter is in the Kondo regime. To leading order in j T , the poles of the 
expression (|67[) are (see Fig. [9] at large d 2 ) 

,,2 A2 ^2 



A li2 = ±iA/t; 2 + A 2 - 



•\ 2 

^2 



2(« 2 + A 2 ) $ 2 ' 



and the amplitudes read 
A 2 



^1,2 — 



2(w 2 + A 2 ) 



O 



A, = 



A3 = 



A 2 



A 2 
-*2 



A? ^ ' 



O 



U 2 



(68) 



(69) 



The expressions (|67l) - (|69p may now be compared with the corresponding ones of a 
Active single biased spin-boson system with parameters A x and v in the white noise 
limit at scaled temperature $. The part that is symmetric in the bias reads 

A 2 + 2dX + v 2 + tf 2 



r (sb) 



(A)). 



A 3 + 2i)X 2 + (A 2 + v 2 + d 2 )X + A 2 tf 



(70) 
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We see that with the identification -d = r y T = K\/ n t) 2 the expressions ([57)1 and (|7D|) are 
quite similar. Observe, however, that the damping rate of the oscillation is somewhat 
different because of the term 2d\ 2 in eq. (JTOJ) instead of i)X 2 in eq. JST]) [2]. Most 
importantly and interestingly in the correspondence, temperature maps on the inverse 
of it. 

7.2. Linear response limit: spin-boson environment as a structured bosonic bath 

We should expect that, in the weak-coupling limit, the spin-boson environment is 
Gaussian and can be represented by a resonant power spectrum of a bath of bosons. 
The Gaussian approximation of the spin-boson environment is found by matching the 
power spectrum of the coupling of the c-spin to the spin-boson environment |14j (with 
normalization as in eq. ([2|) 

v 2 r°° 

S T (u) = —Re / dt([r z (t)r z (0) + rMr z (t)])o^\ (71) 

with that of a harmonic oscillator bath. In the white-noise limit of an unbiased spin, 
the symmetrized equilibrium correlation function Re {T z (t)T z (0)) coincides with the 
expectation (r z ) t . Thus we obtain 

S» = (2/n)v 2 Re{r z (\ = -iu,)) with (r,(A)) = — X+ / 2 , - A2 . (72) 



A(A + tf 2 ) + A 



The resulting power spectrum is that of a structured bath of bosons with a resonance 
of width $ 2 at frequency w = A 2 , 

_ 2tf 2 v 2 A 2 

5>) - — w^w+w ■ ( } 

Due to the coupling to the spin-boson environment, the spin a performs damped 
oscillation, (a z ) t = cos(A 1 i) e _7d ™'. Upon calculating the decoherence rate in order 
v 2 , i.e., the so-called one-boson-exchange contribution of the effective boson bath, we 
obtain 

7dec A tKi) 2[(A|-A?) 2 +l?iA?] ' l ' 

The analysis is completed by observing that this form emerges also upon calculating 
directly 7 dec from the pole equation D 1 (X) = with the form ([M]) . 

8. Conclusions 

We have studied the dynamics of a spin or qubit a coupled to another spin, which 
could be, for instance, another qubit, or a bistable impurity, or a measuring device. 
We have solved the dynamics exactly for white-noise reservoir couplings, and we have 
studied the rich behaviors of the dynamics in diverse limits ranging from weak coupling 
and/or low temperatures to strong coupling and/or high temperature. We have also 
analyzed the effects of a spin-boson environment on the spin dynamics in the Gaussian 
and non-Gaussian domains. 

This paper has not attempted to perform applications to already available 
experiments, instead we have tried to make some general points on complementary 
regimes and on the crossovers in between. 
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One possible simple generalization beyond the white noise limit would be to 
replace the white-noise bath correlations in time intervals in which the Laplace variable 
is irrelevant by the full quantum noise correlation, for instance 

$ f°° 

A 2 2 c -» A 2 c cos(nK c ) ds cos(vs) e~ Q '<^ . (75) 

The advantage would be two-fold: (i) the noise integral is known in analytic form for 
the Ohmic correlation function fbl]) [2|, and (ii) with this substitution the algebraic 
form of pole equation and residua would be left unchanged. 

One further generalization of the Hamiltonian (fTJ) is an applied bias acting on one 
or both of the spins, e.g., of the form e 1 a z and e 2 r z . This important extra ingredient 
could be taken into account exactly in the white noise regime, and would lead to 
additional shifts of the Laplace variable in all blip states of the a- and r-spin. Then 
one would end up with expressions of the form (|32[) and (|47[) with polynomials ramped 
up by bias terms. This extension will be discussed elsewhere. 

Finally, extension of the analysis of the dynamics to the regime T ^ f2 ± requires to 
revert to the original expression ([2|) and compute its effect perturbatively in the one- 
boson-exchange approximation. This can be done either with the self-energy method 
presented in Ref. [5] or with the Redfield approach. One then find, e.g., that the 
actual equilibrium state of (c z t z ) for if, <^C 1 is 



z z / eq 



O 2 



' manh(/3f2/2) - s tanh(/3(5/2) 1 . (76) 



This reduces for T > Cl ± to the previous form ([5^]) found in the white-noise regime. 
The corresponding extension of the analysis and of the results given in Subsection 16.11 
to the domain T < f2 ± will be reported elsewhere. 
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